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Abstract 

We extend the study of Kapitza-Dirac diffraction to the case of two- 
particle systems. Due to the exchange effects the shape and visibility of 
the two-particle detection patterns show important differences for iden- 
tical and distinguishable particles. We also identify a novel quantum 
statistics effect present in momentum space for some values of the initial 
particle momenta, which is associated with different numbers of photon 
absorptions compatible with the final momenta. 

PACS: 42.50.Xa; 03.75.Dg; 61.05.J- 



1 Introduction 

A long time ago Kapitza and Dirac proposed that a standing-wave light field 
can act as a diffraction grating [1] . This proposal is interesting in two main as- 
pects. On the one hand, from the fundamental point of view, it provides a nice 
demonstration of the wave-particle duality where the diffraction grating is not 
massive, but made of massless entities, the photons. On the other hand, from a 
more practical perspective, it allows for the design of matter-wave interferom- 
eters (see, for instance, [2]). The effect has been experimentally observed, first 
with atomic beams [3, 4], later with cold atoms [5], and finally with electrons 
[6, 7]. 

More recently, following the seminal work of Hong-Ou-Mandel (HOM) [8], 
the effects associated with the quantum statistics of the particles have been 
extensively studied in many-particle interferometry. These studies were mainly 
concerned with two-photon interferences originated by the interaction in a beam 
splitter. Later, it has been suggested that the exchange effects can also play 
an important role in interferometry of identical massive two-particle systems 
by diffraction gratings [9]. In particular, in that paper it was shown that the 
diffraction patterns originated at a single slit are very different for distinguish- 
able particles and fcrmions and bosons. 

It seems natural to study the behavior of two-particle systems interacting 
with the Kapitza-Dirac arrangement. However, this is an almost unexplored 
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subject. Up to our knowledge, only in [10] has a numerical simulation of some 
basic aspects of the problem for identical particles been considered. We want to 
present in this paper a more general treatment of the subject for two-particle 
systems, emphasizing the differences with distinguishable particles. We shall 
mainly analyze two aspects of the problem: the spatial dependence of the de- 
tection patterns, and the exchange effects present in momentum space. With 
respect to the first point we shall find that, as in the case of a single slit re- 
ported in [9], the patterns show notorious differences for distinguishable and 
identical particles, in particular, for the shape and visibility of the interference 
figures. A similar behavior is also found for the correlation functions. The sec- 
ond aspect, exchange effects in momentum space, is a less studied subject. We 
shall show that for some values of the initial momenta of the particles a novel 
effect is present in our system. It is reflected in the changes experienced by the 
probabilities of finding the particles in some final states. The physical mecha- 
nism underlying these changes is the possibility of different numbers of photon 
absorptions (the mechanism in the basis of the diffraction of the particles) com- 
patible with the final state of the particles. This effect occurs in addition to the 
usual bunching and antibunching effects in momentum space, which take place 
for close values of the momenta. 

The calculations will be first carried out for single-mode states to present 
the main ideas in a simple way. Later, we shall move to the more realistic case 
of multi-mode states. 

The plan of the paper is as follows. In Sect. 2 we present the basic equations. 
The single-mode diffraction patterns are evaluated in Sect. 3, where we also 
briefly discuss the behavior of correlation functions. Section 4 deals with the 
same problem, but for multi-mode states. In Sect. 5 we move to momentum 
space, in which we describe the existence of novel exchange effects for some 
values of the parameters. Finally, in Sect. 6, we discuss the principal results of 
the paper. 

2 General expressions 

We consider a pair of particles, distinguishable or not, interacting with an op- 
tical standing wave, usually a laser beam. The wave acts as a Kapitza-Dirac 
diffraction grating (see Fig. 1). After the grating we place detectors measuring 
the interference pattern, which can be obtained after many repetitions of the 
experiment. 

The particles passing through the standing-wave light experience a potential 
of the form V = V cos 2 k L x, with fc^ the wavenumber of the light field and x the 
coordinate in the direction parallel to the light beam [11]. The wavefunction 
after the interaction can easily be calculated using standard techniques (see 
[11] for the diffraction and Bragg regimes and [12] for the intermediate region 
between them). For the short interaction times between the particles and the 
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light field the Raman-Nath approximation, equivalent to neglect the particle 
motion during the interaction, holds [2]. This approximation is usual in the 
diffraction regime, the only one we shall consider in this paper. The single-mode 
wavefunction of a particle prior to the interaction is given by ip(x,X,t = 0) = 
exp(i(k x + K X)), with X the coordinate of the direction perpendicular to the 
optical wave, and fco and K the initial wavenumbers. The evolved wavefunction 
after the interaction is given by 

ip(x,X,t) = e- lVt/h i;(x,X,t = 0) = e ^ x e ^v t cos 2 k L x/h^k oX ^ 

where t is the interaction time. The evolution in the X-axis is not modified by 
the interaction. 




Figure 1: Schematic representation of the arrangement. The continuous and 
dashed lines represent the two particles 

Recalling the well-known expressions exp(zzcos^) = J2 n i n J n {z) exp(inip), 
with J„ the Bessel functions, and cos 2 <p = \ + \ cos 2ip, we easily obtain 

oo 

il>(x,X,t) = e iK ° x J2 b n e i{ - 2nkL+k ^ x (2) 

n— — oo 

with b n = i n e- tw J n (-w) and w = V Q t/2H. 

The wavenumber of the particle can only be modified by double recoils. 
For atoms, the first one is associated with the photon absorption whereas the 
second one corresponds to stimulated emission. In the case of electrons, the 
double scattering can be understood as a stimulated Compton scattering [11]. 

We consider now the case of two particles. We denote by ^k K a { x ^ X, t) and 
tp qo Qo(y,Y,t) two wavefunctions, with an obvious notation. When the particles 
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are identical the usual product wavefunction ^(x, X, y, Y, t) = ipk a K ( x i X, t)tp qo Q (y, Y, t), 
valid for distinguishable particles, must be replaced by 

V(x, X, y, Y,t) = -j= (ip ko K (x, X, t)ip qoQo (y, Y, t) ± ip kaKo (y, Y, t)ip qoQo (x, X, t)) 

(3) 

In the double sign expressions the upper one holds for bosons and the lower one 
for fermions. 

From the above expressions one can see that, in the case of identical Particles, 
we have simultaneously two different interference effects. On the one hand, for 
both distinguishable and identical particles, the distributions \ip ko K (x, X, t)\ 2 
display the interference effects associated with the diffraction grating. On the 
other hand, for identical particles we have another interference effect that is 
not present for distinguishable ones. This follows immediately from the term 
±2Re(^l oKo (x 7 X 7 t)ip* qoQo (y,Y,t)tp koKo (y 7 Y,t)ip qoQo (x 7 X 7 t)), contained in the 
expression of \^(x, X, y, Y, t)\ 2 . This form agrees with the standard interpreta- 
tion of exchange effects as interference effects [13]. 



3 Spatial probability distributions 

In this section we evaluate the spatial distribution of simultaneous two-particle 
detections and the correlation functions associated with it. The simplest ex- 
perimental implementation of the arrangement consists of two detectors, one 
fixed at a given position and the other placed at different points in successive 
repetitions of the experiment. In a first step we restrict our calculations to 
single-mode states, postponing the discussion of multi-mode ones to the next 
section. 

The evaluation of the probability distribution is simple. We assume the 
detection time is fixed at t and we can drop the temporal variable from all the 
expressions. We rewrite Eq. (2) as 

ip(x,X) = e iK ° x e ik ° x (t>(x) (4) 

with 

oo 

0(x)= ]T b n e i2nk ^ x (5) 

71— — OO 

With this notation the two-particle probability distribution becomes 

|* dls (x,X,y,y)| 2 = |0(x)| 2 |0(y)| 2 (6) 
for distinguishable particles, and 

|*(x,X,y,F)| 2 - |,/)(x)| 2 |0(2/)| 2 (l±cos((^o-Qo)(^-^)+(fco-?o)(x- 2 /))) (7) 
for identical ones. 
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We assume from now on that X = Y, i. e., the positions of the detectors in 
the direction perpendicular to the light field to be equal. Several consequences 
easily emerge from the above expressions. In the case of identical particles, 
the distribution is the product of the distinguishable distribution by the term 
1 ± cos((fc — qo)(x — y)). The first one reflects the two-particle interferences 
associated with the dispersion by the light field. It is a function, through </>, of 
and b n (V t) the parameters of the optical diffraction grating. That dependence 
is similar for both identical and distinguishable particles. On the other hand, 
the term containing the cosine function is related to the exchange effects. For 
fixed x and y it is only function of the initial momenta of the particles in 
the direction of the light field. The bunching and antibunching effects directly 
emerge from the above equations. In the case of bosons, for x ~ y we have 
I^b^, A, y, A)| 2 s=y 2|</>(x)| 2 |</>(y)| 2 , i. c., the probability of two-boson detection 
almost doubles that of two distinguishable particles. For double detection at the 
same point, x = y, we have, as discussed in [9], |\&b(x, X, x, X)\ 2 = 2|0(x)| 4 , 
i. e., a dependence on the fourth power of the wavefunction modulus. By 
contrast, for fermions we have for x ~ y, | x I'_f(x, X, y, X)\ 2 w that is the 
antibunching effect. When x — y the probability is strictly null according to 
the exclusion principle. There is another situation in which the two-fcrmion 
detection probability is, for any x and y, identically null. It occurs for equal 
initial momenta in the direction of the light field, qo = k$. This is so even when 
Kq 7^ Qo (in the case of Kq = Qo it is evident that the distribution must be 
null because it is impossible to prepare two fermions in the same state). 

Next we present a graphical representation of the above equation. As dis- 
cussed before we fix one of the detectors at y = and move the other at different 
positions x. A simple calculation gives the explicit form for (f>(x) and its squared 
modulus: 

m>n 

\cf>(x)\ 2 = l + 2(-l)" +m J„HJ m Hcos((m-n)(2fc i2; + |)) (8) 

In the calculation we have used the properties J2 n \bn\ 2 = 1 that derives 
directly from the normalization condition, and J n (—w) = (— I)™ J n (w). 

The normalization of ^dis is automatically guaranteed by the normalization 
of <j). On the other hand, for identical particles l^l 2 contains terms of the form 
cos((m — n)2fciX + 7r/2) cos((fco — qa)x). The integration over x of these terms is 
not null (see, for instance, [9] where similar integrations are carried out) when 
(m— n)2kL = ±(fc — qo)- When this condition is fulfilled the normalization is not 
given by |0(O)| 2 £ n \b n \ 2 = I, but by the expression |0(O)| 2 (£„ |6„| 2 ± /) = 1 
with / the result of the above integration. 
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Figure 2: In the vertical axis we represent the two-particle detection probability 
and in the horizontal one the coordinate x. The red, black and blue curves 
correspond, respectively, to bosons, distinguishable particles and fermions. We 
use the values w = 0.2, k = 0.9, qo = —0.9, and fcj, = 1. 

The figure shows a large increase of the visibility of the two-particle detection 
probabilities for identical particles. As usual the visibility is defined as 

v = m^x,o,x)\l ax -Mx,x,o,x)\ 2 min 

In the case of distinguishable particles the probability approximately oscillates 
between 0.98 and 1.02, whereas for identical ones it does between and 2. Then 
we have, respectively, Vdis ~ 0.02 and Vide ~ 1, i. e., a very large difference. 

The curves for bosons and fermions are very similar, with small changes of in- 
tensity due to the modulation introduced by the multiplicative factor |(/>(a;)| 2 |0(O)| 2 
and a phase- n displacement associated with the sign ± (which explains the dif- 
ference between bunching and antibunching at x ~ 0). 

Now, we briefly analyze the behavior of the correlation functions. These 
functions have been extensively used in the experimental study of the (anti)bunching 
effect of free particles released from optical lattices or magnetic traps [14, 15, 16]. 
In particular, in the last two references the existence of periodic correlations 
was shown between density fluctuations of atoms released from optical lattices. 
These correlations reflect the underlying structure of the lattice. 

The correlation function is defined as 

C h) = \j \*( x > X >* + V,X)\ 2 dx (10) 
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with d — 2ir/kL representing the spatial periodicity of the optical diffraction 
grating. In our case, this expression transforms into 

C{rt) = - d {l ± cos(( go - ko)v)) J" |0(^)| 2 |^(^ + v )\ 2 dx (11) 

The evaluation of the integral is simple but lengthy. We only present the final 
result (a similar calculation can be found in [9]): 

C{rf) = (1 ± cos((<7o - ko)rj)) x 

(1 + 2j2(-l) n+m+r+S JnMJ m (w)J r (w)J s (w) cos(2(r - s)k L rj)) (12) 
N 

with Af denoting the set of n, m, r, s with the constrains m > n, s > r and 
m — n = ±(r — s). 

This correlation function is the product of two periodic functions. The first 
one, that in the left side, is only relevant for identical particles. For distinguish- 
able ones it reduces to 1 and, consequently does not vary for different initial 
preparations of the pairs. For identical particles it only depends on their initial 
momenta. For bosons it reaches maximum values for any r] — 2NTr/(q — k ) 
with N an integer and minimum ones for r\ = (27V + l)ir/(q — fc ). A similar 
behavior holds for fermions when an additional phase- it is taken into account. 
On the other hand, the second multiplicative factor is similar for both distin- 
guishable and identical particles and bosons and fermions. It only depends on 
the parameters of the optical grating, w and As all the terms of the type 
(r — s)k[ J are contained in this factor, we have in the correlation function all the 
periods generated by the optical grating. As discussed in [15, 16, 9], it reflects 
the underlying structure of the diffraction grating. 



4 Multi-mode states 

Up to now, we have only considered single-mode states. Now, we move to 
the more realistic case of multi-mode ones. The simplest way to study them 
is to assume that the distribution of initial wavenumbers of each particle is a 
Gaussian one. The one-dimensional Gaussian distribution is given by f(k ) = 
(47r) 1//4 <7~ 1/2 cxp(-(fc - A) 2 /2cr 2 ) with a the width of the distribution and A its 
central value. The initial wavefunction reads as ip(x,t = 0) = J dk f(k )e lk " x . 
Note that, for the sake of simplicity, we have not included the variable A, which 
is irrelevant in the following discussion. The wavefunction after passing through 
the light grating is: 

oo „ 
^(x)= b n e ankLX dk f(k Q y koX 

71 — — OO 

OO 

~e-^ 2 / 2 J2 Kj 2{nkL+K)x = e" lV/ VA(i) (13) 

n— — oo 
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where we have carried out a trivial integration over ko, and i/'A represents the 
wavefunction of a particle with initial wavenumber A. For the matter of simplic- 
ity we have not included the constant coefficients of the distribution and those 
derived from the integration. We conclude that the multi-mode wavefunction 
after the interaction equals that of a single-mode particle with the central value 
of the distribution, but spatially modulated by a Gaussian distribution. 

In the next step we consider two particles in multi-mode states, with central 
values A and T and widths a and ji. The probability distribution after the 
interaction is 

l*(*,y)| 2 = l -e-^e~y 2 ^\4>{x)\ 2 \<t>(y)\ 2 + 

\e^e-^\m?\^)?± 
e- ( * W)(CT2+ ^ )/2 |^)| 2 |<X2/)| 2 cos((x - y)(A - T)) (14) 

We represent this distribution. We take the same values of Fig. 2, in particular, 
A = fco and T = qo. For the width of the distribution we take a 2 = jj 2 = 0.2. 



Figure 3: The same as in Fig. 2, but for multi-mode states. 

As in the single-mode case, there is a notorious difference between the curves 
of distinguishable and identical particles. For the first one, there is an almost flat 
distribution in the center followed by an exponential-like decreasing superposed 
with small oscillations (associated with |0(O)| 2 |0(ir)| 2 ). On the other hand, for 
identical particles there is an interference pattern. In the proximities of the 
point x — we observe the bunching and antibunching effects. The values of 
the visibility are clearly different for the three figures, but not so much as in 
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the single-mode case. These contrasts in the visibility are strongly enhanced for 
small values of a. 

5 Momentum space 

In the two previous sections we have studied the spatial two-particle interference 
patterns. Now we consider the momentum space to see how the exchange effects 
manifest in this representation. The methodology is similar to that used above, 
first we consider single-mode states where the principal ideas can be analyzed 
in a simple way, and later we address the multi-mode case. 

The wavefunction in momentum space (the momentum and wavevnumber 
are related by the trivial relation p = hk, and it is justified to speak about the 
momentum space although really we are dealing with wavevenumbers) is given 
by the Fourier transform of the wavefunction in the position representation, 
$(fc) = J dxe~ lkx Tp{x). For the identical two-particle system we have 

$(k,q) = fdxf dye-^+™^{ Xl y) = ±JS> ko {k)% {q) ± -±=<t> ko (q)<t> qo (k) 

(15) 

where the single-particle wavefunctions are characterized by the value of the 
initial wavenumber (for multi-mode states by the central value of the distribu- 
tion). 

The experimental variable to measure is the probability of detecting one 
particle with the value k and the other with q. In the arrangement we must 
replace the detectors of the previous sections by momentum measurement de- 
vices. In a large series of measurements, for instance fixing k and scanning for 
different values of q, we can compare the experimental data with the theoretical 
distribution: 

i*(mi 2 = |i**o(fc)i 2 i**,(9)r + ^i<M<z)i 2 i<iv(fc)i 2 

±Re(*t (k)*' qo (q)* ko (q)* go W) (16) 

This expression shows that the probability distributions in momentum space 
have a structure similar to that in the position representation. In particular, 
for k s=s q we have the bunching and antibunching effects in momentum space. 
This effect has been numerically discussed in [10]. 

To get a better understanding of that distribution we shall consider specific 
examples, starting with the single-mode states. 

5.1 Single-mode states 

We consider successively the cases of a single particle, two distinguishable Par- 
ticles, and two identical ones. 
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One particle. This case is particularly simple because the wavefunction 
reduces to a superposition of Dirac's deltas: 

$(fc) = b " 6 ( k - 2nk L - M (17) 

n 

For a single particle the probability of finding in a momentum measurement the 
value 2nk^ + fco is \b n \ 2 because, in the sense of the distributions, for r ^ n the 
terms 8(k — 2nk^ — ko)S(k — 2rk^ — fco) are zero. 

Two distinguishable particles. Now, the probability distribution is 

l*fc (fc)| 2 |$ 90 (g)| 2 = \ b M 2 5(k - 2nk L - k )S(q - 2mk L - q Q ) (18) 

The probability of obtaining in simultaneous measurements the values 2nk^+ 
fco and 2mfci + q$, denoted as P(n,m), is | 

To illustrate the form of this probability we graphically Represent it in Fig. 
4 it for the lower values of n and m. Note that this probability is equal to the 
probability of n absorptions by the particle with initial wavenumber fco and m 
by that with qo. For small values of w there is only one absorption. When w 
increases the probability of two absorptions becomes dominant. The probability 
of symmetric absorption (one photon each particle, P(l, 1)) is much larger than 
the asymmetric one (one particle two photons and the other none, P(0, 2)). For 
larger values of w the terms with three absorptions become important. However, 
here, there is again a notorious difference between P(l,2) that increases with 
w and P(0, 3) that remains negligible for all the range of values considered. At 
the end of the graphic, the term P(2, 2) becomes comparable to (or larger than) 
the other terms. 
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Figure 4: Representation of the probability P(n,m) versus w. The brown curves 
correspond to P(0, 1) (continuous), P(0, 2) (dashed), and P(0, 3) (dotted), the 
green ones to P(l, 1) (continuous), and P(l, 2) (dashed), and the purple one to 
P(2,2). 

Two identical particles. The two direct terms in |$(fc,<7)| 2 are similar to 
those for distinguishable particles. We evaluate now the term associated with 
the exchange effects: 

fa)**. (?)**(*) = 

b* n b* m b r b s 5(k - 2nk L - k ) x 

n,m,r,s 

6(q - 2mk L - q )5(q - 2rk L - k Q )S(k - 2sk L - q Q ) (19) 

The product of the two deltas containing k is zero (in the sense of the distribu- 
tions) unless the relation 2nkj^ + fco — 2skj^ + qo holds. Similarly, we must have 
2mkh + qo = 2rk^ + fco hi order for the product of the deltas containing q not 
be null. These two relations can be expressed as 

n- s = N ; m-r - -N ; N = ^— — (20) 

2k L y ' 

with N an integer. The difference between the initial wavenumbers must be an 
integer number of times 2/cl. 

When these relations hold, we have that the terms related to the deltas 
containing n and m must be rewritten as 

\\b n bm\ 2 S{k - 2nk L ~ k )S(q ~ 2mk L - q ) + 
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-^\b n b m \ 2 S(q - 2nk L - k )S(k - 2mk L - q ) ± 
Re(b* n b* m b m+N b n - N )5(k - 2nk L - k )S(q - 2mk L - q ) (21) 

The two first terms correspond to the absorption of n photons by a particle with 
initial wavenumber kg and m by one with qo. If (<Zo — fco)/2fc.L is not an integer, 
the third term does not contribute and we have the two final wavenumbers 
2nk^ + fco and 2mfci + qo with probability |6„& m | 2 . On the other hand, when 
the condition holds we obtain the same final wavevenumbers 2nk^ + fco and 
2mfci + qo, but with probability 

P N (n, m) = \b n b m \ 2 ± Re{b* n b* m b m+N b n - N ) (22) 

This result can easily be understood in terms of indistinguishability of alterna- 
tives. If the final wavenumbers are the same for both types of absorptions, we 
cannot distinguish if the final result corresponds to the alternative of absorp- 
tions n and m, or to the other alternative, n — N and m + N. In quantum theory 
the amplitudes of probability for indistinguishable alternatives must be added, 
obtaining an interference effect. The indistinguishable alternatives correspond 
to different numbers of absorptions yielding the same final wavenumbers. The 
situation is different for distinguishable particles. In this case, we can know in 
principle at any instant if the particle is of one type (that whose wavenumber 
is denoted by k) or the other (wavenumber q) and, consequently, if its initial 
wavenumber was k or q . Then one particle of the first kind cannot reach 
the final wavenumber 2nkh + fco by n — N absorptions. There are not differ- 
ent alternatives to reach the final wavenumber and, there is not an interference 
effect. 
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Figure 5: Representation of the exchange effects. Black and white circles rep- 
resent particles labeled by k and q. (a) The standard (anti)bunching effect in 
the spatial representation. The detections can occur in two alternative ways, 
these represented by the continuous and dashed lines, (b) The same effect in 
momentum space. The dashed horizontal lines correspond to the two momenta 
at which the particles are detected and the continuous vertical ones to the tra- 
jectories in momentum space. The (anti)bunching is observed when the two 
final momenta are very close, (c) The effect discussed in this paper. The dotted 
lines correspond to the initial momenta. The arrows represent the absorption 
of photons. Being the particles identical the two represented alternatives are 
indistinguishable . 



We remark that this novel exchange effect is different from (anti)bunching 
(see Fig. 5). The last one only takes place when we consider particles very close 
in momentum space. In contrast, the effect described here only depends on the 
initial values of the momenta (which can be rather different for all the cases 
with N 7^ 0) and can occur for well-separated final momenta. 

The previously described effect is a purely quantum one. To justify this 
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statement we show that it cannot be described by a classical treatment. We use a 
qualitative model where the particles collide with the photons and the collisions 
are ruled by the classical law of momentum conservation (classical description 
of the Compton effect assuming the existence of photons). We assume that 
all of the momentum of the photon is transferred to the particle, in such a 
way that its momentum changes from p to po + hki,. In the cases with 2n 
collisions the final momentum of the particle would be the same predicted by 
quantum theory. However, in the absence of additional unnatural assumptions, 
the classical theory allows for odd numbers of collisions. The agreement between 
the model and quantum theory would be even worse for couples of particles. 
Classically, both particles behave in an independent way, and the stochastic 
collision events are uncorrelated. The probability of one particle experiencing 
2n collisions and the other 2m factorizes. The result would be similar to that 
for quantum distinguishable particles. In conclusion, in the classical model (a 
reasonable one) there is not room for the previous exchange effect. 

Next we represent these probabilities for a particular example. A simple 
calculation using the expression for b n and the property of the Bessel functions 
J_ n (u>) = (~l) n J n (w) gives P (1,0) = J2J 2 ± Jpl P x (1,0) = J2J 2 ± j\Jl 
and P_i(l,0) = JfJ§ T J\ J Q J i- We see that P (1,0) = Pi (1,0). 




0.0 0.5 1.0 1.5 



Figure 6: Representation of the probability P of finding the final wavenumbers 
2nkL+ko and 2mkL+qo versus w. The black curve represents the distinguishable 
case for n = 1 and m — 0. The red and blue ones correspond to bosons and 
fermions with N = 1 (continuous) and N = — 1 (dashed). 

The figure clearly shows that the probabilities of finding the particle in the 
targeted final wavenumber state are different for the all the cases. For N = 1 
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the detection rate notoriously increases for bosons, whereas it becomes null for 
fermions. The detection rate is not always enhanced for bosons; for N = — 1 it 
is smaller than that of distinguishable particles. Similarly, the fermion rate can 
increase with respect to that of distinguishable particles (also N = — 1). 



5.2 Multi-mode states 

In this subsection we analyze how the above results are modified in the more 
realistic case of multi-mode distributions. We mainly restrict our considerations, 
as in Sect. 4, to Gaussian distributions which can be tackled analytically. As in 
the previous subsection we consider successively the cases of a single particle, 
two distinguishable particles, and two identical ones. 

One particle. The wavefunction in momentum space for a particle with a 
multi-mode distribution f(k ) of initial wavenumbers is easily obtained: 

I dk f(k )b n S(k - 2nk L - k ) = ^ b n f(k - 2nk L ) (23) 

n n 

The probability of finding the particle with final wavenumber k is |$(fc)| 2 = 
S n m b%bmf*{k — 2nkL,)f{k — 2mkL)- If the mode distribution is a Dirac's delta 
(the single-mode distribution), /(fc ) = ^(^o ~ K ), the probability becomes 
\3>{k)\ 2 =Y, n \ h n\ 2 8{k-2nk L -K ). 

When the initial distribution is a Gaussian, /(fco) ~ e -( fc o-A) /2a ^ we nave 

|$(fc)| 2 ~ ^me" (fc " M '" A ' : l /2,\-(k-2 mkL -Kf '/2a 2 (24) 

These two exponentials correspond to two Gaussians with the same width of 
the initial ones, but centered around 2nk L + A and 2mk L + A. Two different 
regimes can be obtained. For 2k^ 3> cr we have that the overlapping between 
the Gaussians centered around 2nkh + A and 2{n ± l)k^ + A (and, of course, 
any m with \n ± m\ > 1) is negligible. The sum in the above expression re- 
duces to the diagonal terms, |$(fc)| 2 ~ J2 n \b n \ 2 er {k - 2nk t- K f I* 2 . After the 
interaction, we can only detect the particle with wavenumbers contained in the 
Gaussian distributions centered in the points 2nk^ + A (with \j\[2 times the 
width of the initial one). There is not a contribution of the crossed terms to the 
probability of detection. In the limit of very peaked Gaussians we recover the 
behavior described by Dirac's delta distributions. The other regime takes place 
when a > 2kL- In this case, there is a non-negligible overlapping between the 
distributions centered around 2nkL + A and 2(n ± l)fcc + A (and, perhaps, the 
other m with \n — m\ > 1). The crossed terms can no longer be neglected, lead- 
ing to interference terms. These interference terms can be understood by the 
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impossibility of distinguishing if a particle detected with wavenumber k belongs 
to one or the other of the distributions. 

Two distinguishable particles. The discussion closely follows that of 
one particle. If the initial distributions of the particles are /(fco) and <?(&o) the 
probability of finding one particle with k and the other with q is 

|%(fc)| 2 |$ 9 (g)| 2 ~ Yl b* n b m b* r b s f*(k-2nk L )f(k-2mk L )g*(q-2rk L )g(q-2sk L ) 

(25) 

We assume the two distributions to be Gaussian ones. For 2k l ~> <? and 2k^ >• 
\x we have that the overlappings between the Gaussians centered around 2nk^+A 
and 2(n±l)fcL+A and 2nk^ + T and 2(n±l)fc£ + T are negligible, becoming the 
probability ^(fc)! 2 !^)! 2 ~ £„ m \b n b m \ 2 e^ k - 2nk --^ 2 1" 2 e-^- 2nk ^ 2 . 
The probability is the product of the one-particle distributions without crossed 
terms. When the conditions a > 2k^ and/or \i > 2k^ hold we have crossed 
terms leading to interference effects. 

Two identical particles. The two-particle probability contains the terms 
of the form |$/(fc)| 2 |$ 9 (q)| 2 , which can be treated as before, and the exchange 
term that transforms into Ylm m r s R e (K.f*(k ~ 2nkL)b* m g* (q - 2mk L )b r f(q - 
2rkL)b s g(k — 2sfci)). We again consider Gaussian distributions, for which the 
exchange term reads proportional to 

£ Re(b* n b* m b r b s )e-^ k ^ 2 ^ 2 x 

e -(q-2rnk L -Tf /2^ 2 e -(q-2rk L -A) 2 /2a 2 e -(k-2sk L -T) 2 /2^ 2 ^6) 

For the sake of simplicity we take u = (i. For \2(n — s)kL + A — T)| » a or 
1 2 (to — r)fcL — A + T)| » g (for any n and s and to and r) the overlapping 
between the curves is negligible and the product of the two distributions with 
the same argument (k or q) is almost null. The contribution of the exchange 
terms can be neglected. In contrast, when \2(n — s)kL + A — T)| < a for some 
n and s and |2(m — r)k^ — A + T)| < a for some m and r the product of the 
distributions cannot be neglected. The exchange effects become important in 
this case. However, the conditions for the presence of the exchange effect are 
much less stringent than in the case of single-mode states. 

6 Discussion 

We have analyzed in this work the extension of the Kapitza-Dirac effect to the 
case of two-particle systems. The spatial two-particle detection patterns display 
notorious differences for distinguishable particles and fermions and bosons, in 
particular for the shape and visibility of the interference figure. 

We have also shown the existence of a novel exchange effect in momentum 
space. The effect, which modifies the distribution of particles detected with 
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given momenta, only occurs for some values of the initial momenta of the parti- 
cles (or in the case of multi-mode states for some values of the parameters of the 
initial momenta distributions) . For multi-mode states these conditions are much 
less stringent than for single-mode ones. The verification of this effect would be 
interesting in several aspects. Exchanges effects are a striking manifestation of 
the departure between classical and quantum descriptions of physical systems. 
Any new example of these differences is worth investigating. The novel effect 
has no relation with other exchange effects such as (anti)bunching, exclusion in 
atoms, or degeneracy in gases. The physical underlying mechanisms are differ- 
ent, in our case, different numbers of photon absorptions. Moreover, the confir- 
mation of the effect would corroborate the validity of the (anti)symmetrization 
principle in a framework where strong interactions with other types of particles 
are present. 

In this paper, we have only considered the spatial part of the wavefunction. 
This is equivalent to assume that the spin states (and the electronic states for 
atoms) are symmetric and do not play any role in the problem. The extension 
to the case of antisymmetric spin or electronic states, where the relative sign of 
the spatial part of the wavefunction can be reversed, is simple. 

Some other aspects of our approach must be commented on. The wavefunc- 
tions have just been evaluated after the interaction with the optical grating. 
The detectors must be placed at these positions, and consequently, we are in 
the near-field regime. However, in this type of problems one usually works in the 
far-field one [7] . In the free evolution of the particles between both regimes there 
is a dispersion of the wavefunctions. Fortunately, that evolution is simple (free 
evolution) and can easily be described. For instance, in the spatial picture we 
have that the width of the Gaussian packet increases with time. As the momen- 
tum peaks are determined via spatial detection [7] we must take into account 
that broadening to correctly interpret the data. For well-separated peaks the ef- 
fects of the spreading are negligible. For peaks with appreciable overlapping the 
range of values of the initial momenta for which the exchange effect described in 
this paper takes place increases (as we can easily see with an argument similar 
to that used for multi-mode states in subsection 5.2). An exact treatment of the 
problem would require a quantitative evaluation of the evolution of the multi- 
mode states, instead of the qualitative one presented here. However, this will be 
presented elsewhere. Another important aspect of the problem is the question 
of the coherence between the two particles. As is well-known, when they are 
only partially coherent the two-particle interference properties, in particular the 
visibility, are modified. Thus, the relative coherence of the two particles should 
be checked for every type of source of pairs of particles. When it is only partial 
we should introduce the necessary modifications in the formalism. In our pro- 
posal, there is an additional question similar to that of the partial coherence: 
we must analyze if the overlapping between the wavefunctions is complete or 
partial (see also the next paragraph). In the first case we can use a completely 
(anti)symmetrized wavefunction. In the second one, we should include in the 
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wavefunction (or in a density matrix) the property of partial overlapping. This 
is a subject worth analyzing since, as far as we know, it is not considered in the 
literature. Here, by a matter of simplicity, we have assumed complete coherence 
between both particles in the two senses, coherence and overlapping. 

The last question we briefly address here is the possibility of experimentally 
testing the above effects. In the space representation, the interchange effects are 
present when there is a non-negligible overlapping between the wavefunctions of 
the two particles. Then we must prepare the particles in such a way that they 
have a non-negligible overlapping at the time they reach the optical grating. 
Several sources of identical particles, such as Bose-Einstein condensates, optical 
lattices, or magnetic traps have recently been used to study correlation functions 
[14, 15, 16]. If we would be able to select the cases in which only two particles 
are released from any of the above devices we would have an efficient source for 
our problem. We would also carefully test that the times of the arrival of the two 
particles to the optical grating are close enough. In terms of wavepacket spread, 
the peaks of the two probability distributions must reach the optical grating 
at the same time. On the other hand, the interaction strength is given, as in 
the one- particle case, by the parameter w. Using the same values of the single 
particle arrangement [11] we can reach an appreciable value for the number of 
pairs diffracted. 
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